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Q , Abstract. In this paper, we review some of the properties of dense molecular cloud cores. The results 

^ ■ presented here rely on three-dimensional numerical simulations of isothermal, magnetized, turbulent, and 

' self-gravitating molecular clouds (MCs) in which dense core form as a consequence of the gravo-turbulent 

, fragmentation of the clouds. In particular we discuss issues related to the mass spectrum of the cores, their 

lifetimes and their virial balance. 
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, 1 The Simulations 

m 
m 

' We performed 3D numerical simulations of magnetized, self-gravitating, and turbulent isothermal MCs using 
^ ! the TVD code (Kim et al. 1999) on grids with 256^ and 512^ cells (Vazquez-Semadeni et al. 2005; Dib et 
' al. 2007a; Dib et al. 2008a). The basic features of these simulations are: Turbulence is driven until it is fully 
. developed (at least for 2 crossing timescales) before gravity is turned on. The Poisson equation is solved to 
account for the self-gravity of the gas using a standard Fourier algorithm. Turbulence is constantly driven in 
the simulation box following the algorithm of Stone et al. (1998). The kinetic energy input rate is adjusted such 
' as to maintain a constant rms sonic Mach number Mg = 10. Kinetic energy is injected at large scales, in the 
wave number range fc = 1 — 2. Periodic boundary conditions are used in the three directions. In physical units, 
the simulations have a linear size of 4 pc, an average number density of n = 500 cm^"^, a temperature of 11.4 K, 
a sound speed of 0.2 km s~^, and an initial rms velocity of 2 km s~^. The Jeans number of the box is Jtox = 4 
(number of Jeans masses in the box is Mtox/Mj^tox — J^ox — 64, where M^ox = 1887 Mq). The simulations vary 
by the strength of the magnetic field in the box with Bq = 0, 4.6, 14.5, and 45.8 /LtG for the non-magnetized, the 
strongly supercritical, the mildly supercritical, and the subcritical cloud models, respectively. Correspondingly, 
the plasma beta and mass-to- magnetic flux (normalized for the critical value for collapse M/cj) = (47r^G)~^/^; 
Nakano & Nakamura 1978) values of the box are fip^box — oo, 1, 0.1, and 0.01, and ^Mhox = oo, 8.8, 2.8, and 0.9, 
respectively. Cores are identified using a clump-finding algorithm that is based on a density threshold criterion 
and a friend-of-friend approach as described in Dib et al. (2007a). We restrict our selection of cores to epochs 
where the Truelove criterion (Truelove et al. 1997) is not violated in any of them. Thus, the derived properties 
of our numerical cores can be best compared to those of starless prestellar cores. 
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Fig. 1. Core mass function for cores identified at density 

thresholds of 20 and 50 times the average density (full 
thick and dot-dashed line respectively). Over-plotted to 
the CMF is a lognormal fit (full line) and power law fits 
in the intermediate to high mass regimes (dashed lines). 



Fig. 2. Slope of the CMF in the intermediate and high 
mass regimes as a function of the density threshold. The 

slope of the CMF is steeper when cores are identified with 
higher density tracers/thresholds. 



One important area of modern astrophysics revolves around characterizing the initial mass function of stars 
(i.e. mass function of stars at their birth), and assessing whether it is universal (i.e. independent of environment) 
or not. Early determinations of the IMF have been obtained by Salpeter (1955) who showed that the field star 
population can be described by a power law of the form dN/dM = with a ~ —2.35 in the intermediate 
to high mass regimes. Subsequent determinations of the IMF, particularly that of stellar clusters, have derived 
values of a that are in the range ~ [—1.8,-2.7] (e.g., Massey et al. 1995; Sharma et al. 2007). A related 
important issue is that of the relationship between the IMF and the core mass function (CMF). Observations of 
dense cores in nearby star forming regions tend to indicate that their mass distribution is not very different from 
that of the IMF (e.g., Motte et al. 1998; Johnstone & Bally 2006) which would indicate that the shape of the IMF 
might be already set in the early gaseous phase. Using our above described simulations, we have constructed 
the CMF at several density thresholds. Fig. 1 displays the CMF for cores identified at the density thresholds 
of nthr = 20n = 10^ cm~'^ and 50n = 2.5 x 10** cm~^. The derived values of the slopes in the intermediate 
to high mass regimes show a steepening of the slope of the CMF with increasing threshold. Equivalently, the 
width of the CMF decreases with increasing threshold when the CMFs arc fitted with a lognormal distribution. 
Our statistics of cores does not allow us to probe the values of the slopes for cores identified at higher density 
thresholds, e.g., in the range of 200n = 10^ cm~^ at which all cores are gravitationally bound. However, an 
extrapolation of the slope-threshold relation in Fig. 2 seems to indicate that values of the slopes in the range 
[-2.5,-3] can be expected for these threshold values. Nonetheless, it shoud be noted that the power law shape 
(or lognormal) of the CMF can be substantially modified as time evolves under the influence of various physical 
processes such as gas accretion and core coalescence (Dib et al. 2007b; Dib 2007; Dib et al. 2008a,b,c) 



3 Virial Balance 

In Dib et al. (2007a) and Dib & Kim (2007), we investigated the detailed virial balance of the cores formed in 
the simulations. By calculating all the terms of the virial theorem, we could show that the cores are dynamical, 
out of equilibrium structures. In each model, there is a mixed populations of gravitationally bound cores, 
cores that are bound by external compressions, and imbotmd dispersing objects. We compare the diagnostic 
of the gravitational boundedness of cores using the detailed virial balance to that made using other simplified 
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Fig. 3. The virial parameter as a function of tlie mass for tlie clumps and cores in tlie mildly magnetically supercritical 
simulation after t = 1.6 Myr of turning on self-gravity. 

indicators, such as the virial parameter a^ir which compares the sum of the kinetic and thermal energy to 
gravity, the mass to magnetic flux ratio which compares the importance of magnetic support against gravity, 
and the Jeans number which compares the importance of thermal support against gravity. In all simulations, we 
find a trend in which, for gravitationally bound cores, the virial balance of such objects indicates that their inner 
parts (i.e., when cores are defined at high density thresholds) are more gravitationally bound than when the 
same objects are identified with lower density tracers (inner -|-outer parts). The simplified indicators seemed to 
indicate a different result in which the extended objects are more gravitationally bound than their inner parts. 
The discrepancy is nonetheless not very surprising considering the fact that the simplified indicators describe 
only a part of the energetic balance of the cores and that all of them neglect the surface energy terms that 
appear in the virial theorem. A comparison of the mass-virial parameter relation = Mj^ for the different 
simulated clouds (and at different epochs) with the observations (e.g. Bertoldi & McKee 1992; Williams et al. 
2000) shows that the closest agreement with the observations is obtained for the case of the mildly magnetically 
supercritical cloud model. The value of /3 in this model is found to be ^ —0.6 (Fig. 3) in contrast to values 
of /? > —0.5 for the strongly supercritical and non-magnetic cases. This seems to indicate that real molecular 
clouds in nearby star forming regions and in which cores are observed to form might be in a state of near 
magnetic criticality or slight super-criticality. 

4 Cores Lifetimes 

In Vazquez- Semadeni et al. (2005) the time evolution of the cores was studied for simulations of different 
initial mass-to-magnetic flux ratio. Galvan-Madrid et al. (2007) performed mesurements of the core lifetimes 
as a function of the density threshold rithr used to define the cores for the likely case of MCs that are slightly 
magnetically supercritical (see above). The prestellar lifetimes of the cores ranged between a few to several 
free-fall times, with a mean value of Tpre ~ 6 t//, where tff = (37r/32G/9)^/^. Fig. 4 displays the lifetimes of the 
cores in the simulations against their threshold density (which is similar to the initial mean volume density). 
The observational data compiled by Ward-Thompson et al. (2007) is overplotted (filled squares). There is good 
agreement between the prestellar lifetimes of the simulations and the observations. Galvan-Madrid et al. (2007) 
also performed indirect estimations of the starless-to-protostcUar core ratio. They found that it decreases from 
~ 5 at rithr ~ 10'* cm~^ to ~ 1 at iithr ~ 10^ cm~^, in rough agreement with observational estimates (e.g., Lee 
& Myers 1999; Hatchell et al. 2007). 

S. D. would like to thank the organizers of the SF2A-PCMI Dynamique workshop for the invitation to speak. 
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Lifetime vs volume density 
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Fig. 4. Duration of the prestellar stage for the collapsing cores in the simulations as a function of their mean density. 
Diamonds, triangles, and open squares correspond to different runs. The loci of Tpre — 1 iff and Tpre ~ 10 t// are 
marked by the lower and upper long-dashed lines respectively. The linear fit to the simulation measurements (rpre — 6 
tff) is marked by the dash-dotted line. Over-plotted to the simulations are observational data points (filled squares, 
taken from Ward- Thompson et al. 2007). 
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